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We present a unified formalism for describing chemical reaction rates of trapped, ultracold 
molecules. This formalism reduces the scattering to its essential features, namely, a propagation of 
the reactant molecules through a gauntlet of long-range forces before they ultimately encounter one 
another, followed by a probability for the reaction to occur once they do. In this way, the electric- 
field dependence should be readily parametrized in terms of a pair of fitting parameters (along with 
a Ce coefficient) for each asymptotic value of partial wave quantum numbers \L,M). From this, 
the electric field dependence of the collision rates follows automatically. We present examples for 
reactive species such as KRb, and non-reactive species, such as RbCs. 



PACS numbers: 34.50.Cx,34.50.Lf 

Ultracold molecules present researchers with unique 
physical systems that arc a curious mixture of small 
and large energies, and of tiny and enormous length 
scales. Thanks to recent experimental advances, certain 
molecules can be prepared in specific hyperfine states, 
even though they are separated in energy only by ~ 10~^ 
eV [il-Q- Yet, upon colliding, these same molecules ex- 
plore the ^ 1 eV energies afforded by their electronic 
potential energy surfaces. Similarly, the translational ki- 
netic energy of these molecules, set by their temperature 
T, can be as small as ksT ^ 10""'^^ eV. At these energies 
the force the molecules exert on one another can be signif- 
icant on length scales that are orders of magnitude larger 
than the molecules themselves. Because long-range dipo- 
lar forces are experimentally controllable, there has been 
much discussion about the prospect of either limiting or 
enhancing chemical reaction rates by the simple artifice 
of changing an electric field Q. Indeed, an effect of this 
kind has been dramatically demonstrated and explained 
recently Q- Collisions of ultracold molecules are thus, 
in principle, remarkably complicated systems to under- 
stand in detail. On the one hand, every degree of free- 
dom is involved, from the hyperfine states in which the 
molecules are prepared to the complete re-arrangement 
of molecules in a chemical reaction. On the other hand, 
the actual number of observables may be rather small, 
consisting perhaps of a handful of loss rate coefficients. 
The long path connecting complex molecular Hamiltoni- 
ans to what is actually seen in the lab may indeed prove 
intractable from ah initio theory. 

For this reason, it is worthwhile to find simple formu- 
lations of collision theory at ultralow temperatures, espe- 
cially formulations that naturally take advantage of the 
vast differences in energy and length scales present [y, 0] ■ 
Recently, two complementary approaches have accounted 
fairly well for experiments that have observed ultracold 



chemical reactions of fermionic KRb molecules. In the 
first, a multichannel quantum defect theory (MQDT) ap- 
proach has successfully replaced the short-range physics 
by suitably parametrized boundary conditions that ac- 
knowledge both a scattering phase shift and the prob- 
ability of chemical reaction [8|. The boundary condi- 
tions were matched to highly accurate solutions of the 
long-range scattering, which in fact were carefully char- 
acterized analytically, allowing for simple analytic for- 
mulas for scattering observables Q. The second, "quan- 
tum threshold" (QT) approach focused on the fact that 
the molecules had to tunnel through a centrifugal barrier 
with a given probability, which varied with energy in ac- 
cord with the Wigner threshold laws [^] . By fioating the 
value of the tunneling probability at the barrier's peak, 
this method was able to describe in an analytic way the 
chemical reaction probability even in the presence of an 
electric field that polarized the molecules Q. 

In this Rapid Communication we merge the ideas be- 
hind these approaches to arrive at a consistent theory 
of ultracold polar molecule collisions. We will exploit 
the short-range parametrization already afforded by the 
zero-field MQDT approach, complemented by a numer- 
ical treatment of long-range wave function propagation. 
One main result is the classification of molecules accord- 
ing to whether their scattering is universal, with loss 
rates that depend only on purely long-range features of 
the potential energy surfaces; or else non-universal, con- 
taining resonances that carry more detailed information 
about specifics of the interactions. These kinds of field- 
dependent resonances were reported previously |10l4l6l |. 
We find that the locations and contrast of these features 
are specified once the MQDT parameters of the short- 
range physics are given. Therefore, a whole swath of the 
electric-field-dependent coUisional spectrum may be sim- 
ply characterized. 



We begin with the Hamiltonian for interaction of two 
dipolar molecules in well-defined single internal states, 
with reduced mass fj. and intermolecular separation R: 



H = Tr + Vs, + yccnt + K 



vdW 



Vm- 



where k is the incident wave number, and g = 1,2 ac- 
cording as the particles are indistinguishable or distin- 
guishable in their initial channel. 

The S-matrix is characterized by a complex phase 
(1) shift via S = e^*'' , which defines the complex, energy- 



Here, Tr is the radial kinetic energy; l^ont = fi^L{L + 
1)/2^R^ is the centrifugal energy corresponding to par- 
tial wave L; V^dw = —Cq/R^ is the van der Waals in- 
teraction between two molecules, here assumed isotropic; 
and Vdd = C^{L,L'] Ml)/R^ is the dipole-dipole interac- 
tion between the molecules, which couples different par- 
tial waves, but which preserves the projection Ml of this 
angular momentum onto the field axis 19[. These three 
terms identify the long-range physics, denoted collec- 
tively as V\-c. In addition, T4r incorporates all short-range 
physics, such as elastic and inelastic scattering, possible 
resonances to ro-vibrationally excited molecular states, 
or even chemical reaction. We will not deal explicitly 
with Vsr in what follows. 

For the present, we are interested in collisions that 
may result in chemical reactions, rather than hyperfine- 
changing collisions. We therefore ignore all other hyper- 
fine states besides the incident ones. (The theory can be 
adapted to include these later, however.) In our model, 
scattering via the Hamiltonian ([1]) is then a multi-channel 
problem, where the channels are defined by the partial 
wave quantum numbers L. Higher values of L generate 
higher centrifugal barriers, and thus inhibit the passage 
of the molecules to short range where they can react. 
Therefore, only a handful of L's are necessary to de- 
scribe chemical quenching phenomena at ultralow tem- 
peratures. In fact, we consider a single potential that is 
constructed by diagonalizing Vy-c in the partial wave basis 
at each value of R, in the spirit of the Born-Oppenheimer 
approximation, as in Ref. [l6[. 

We therefore reduce the problem to scattering in a sin- 
gle potential, albeit one from which wave function flux 
can leak at small values of R. Within this model, the 
scattering consists of three parts: 1) molecules approach 
one another and transmit some fraction of their inci- 
dent fiux through the long-range potential T4r, to arrive 
at an intermediate separation Rq\ 2) the molecules en- 
ter into the near zone, where they may either react (in 
which case fiux is lost) or else scatter back into the chan- 
nel from which they came, generating a phase shift; 3) 
what's left of the molecular fiux proceeds to infinite R 
and counts toward elastic scattering. In either event, 
scattering is defined via the diagonal scattering matrix 
element in our potential, S, whose magnitude may be 
less that unity if a reaction has occurred. Quite gener- 
ally, elastic and quenching scattering rate constants are 
given, respectively, by 



IC\E)=g^-^\l-Sm' 
K^''{E)^g^{l-\S{E)f), 



dependent scattering length 



l{k) =a{k) -ij3{k) = - 



tan77(fc) 



(3) 



The real power of the quantum defect approach is that it 
provides analytic expressions for the c omp lex scattering 
length for zero-electric- field collisions [Ijl. This follows 
from a careful parametrization of standard wave func- 
tions in the long-range potential, which is assumed to 
consist solely of van der Waals plus centrifugal poten- 
tials. As was shown in Ref. [^, the scattering lengths for 
the lowest partial waves simplify, in the limit fca <g; 1, to 



aL=a 



ay- 



l + (l-s)2 



(2^=1 = —2ai(ka) 



i + y{l - s) 
^ y + i{s-l) 
ys + i{s — 2) 



(4) 
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Here several scale parameters are used, such as 
the Gribakin-Flambaum mean scattering length a = 
27r(2/xC6//i'^)^/'*/r(l/4)^, and its p-wave analogue oi = 
ar(l/4)6/(1447r2r(3/4)2) = 1.064a. The parameters 
that are specific to each particular scattering problem 
are instead the real part of the zero-energy scattering 
length, a, also given in its reduced form s = a/a; and 
the effective short-range channel coupling strength y. In 
the quantum-defect point of view, y stands for the prob- 
ability of chemical reaction once the molecules get close 
together: when y = 0, chemical reactions are forbidden, 
and the scattering is purely elastic; whereas when y = 1, 
reactions occur with maximum probability [8|. 

Thus the parametrization of scattering observables fol- 
lows as given above, whereas the actual values of the pa- 
rameters s and y will vary from one molecule to the next, 
and may be determined by fitting experimental data. 
Within the theory, their values follow ultimately from the 
value and derivative of the total wave function ip(Rm) at 
a matching radius Rm- In MQDT, Rm represents the 
boundary between long and short-range, or equivalently 
small and large energy scales in the relative motion of the 
molecules. Its value is conveniently chosen to be smaller 
than the characteristic length a, yet larger than the scale 
of any significant short-range physics or chemistry. 

We come now to the main point of this article. At 
the matching point R„i, the interaction potential is suf- 
ficiently deep that the potential, and more importantly 
the wave function 'ip[R„i), are unaffected by turning the 
electric field on. This refiects the physical fact that lab- 
oratory strength fields, inducing energies on the order of 
10~^ eV by polarizing the molecules, has no effect what- 
soever on the eV-scale chemical reaction processes, or the 
perhaps 0.1 cV depth of the interaction potential at Rm- 
Therefore the parameters s and y can be defined once, at 
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FIG. 1: (Color online) Probability to tunnel to short range, 
versus collision energy, for s- and p-wave collisions of KRb 
molecules, assuming zero electric field and Ce = 16130 a.u. 
[15| . The dashed lines are the analytic, low-energy approxima- 
tions 1 — exp(— 4fc/?). Inset: Dependence of the transmission 
probability, evaluated at the height of the centrifugal barrier, 
as a function of induced dipole moment d. 



zero electric field, and then used at all subsequent, higher 
fields. It only remains to propagate the wave function 
■0 from Rra out to infinity. This part of the process is 
necessarily numerical, since the MQDT parameters for a 
mixed van der Waals-plus-dipole interaction are not yet 
characterized, and in any event are likely to be available 
only numerically. 

In this way the entire field-dependent scattering be- 
havior in contemporary experiments can be succinctly 
summarized. To set this discussion in context, wc first 
consider the transfer function, i.e., the probability that 
the incident flux reaches Rm at a-U, in zero electric 
field. This function is given by the loss probability 
-Pioss = 1 — I'S'P, evaluated for unit short-range absorp- 
tion probability y ~ \. It is shown in Fig. 1 for cold 
collisions of KRb molecules in zero field, for both s- and 
p-wave collisions. In both cases, the low-energy behav- 
ior of Pioss must follow the Wigner threshold laws. In 
the higher-energy limit, P\as& asymptotes to unity, since 
the molecules can then barrel past the comparatively 
weak long-range forces. The dashed lines show approxi- 
mate transfer functions that incorporate the lowest-order 
complex scattering length from Q, namely, /3o = a for 
s-waves and /?i = ai(fca)^ for p- waves. These approx- 
imations arc adequate up to the characteristic energies 
So = h^/2^ia^ for s-wavcs, and Ei = {m^ /{2^^r^CQ)Yl'^ 
for p-wavcs (corresponding to the height of the p-wavc 
centrifugal barrier). For example, for the KRb molecules 
considered in Fig. 1, i?o = 98/iK and Ei = 24.3/iK. 

In the very low energy limit for p-waves, the quan- 
tum threshold model in Ref. [9| also yields the cor- 
rect low-energy behavior of Pioss, and indeed is based 
on the Wigner laws. The only remaining ingredient 
within this model is to normalize Pioss to its correct value 
at the height of the p-wave centrifugal barrier. This 



value, Pioss(Vf,) = 0.37, is also indicated in the figure, 
and it is a universal value that is independent of the 
specific Ce or reduced mass of the collision partners. 
In the QT model the transmission function is given as 
PjQgg = Pyass{Vb){E /ViY^'^ . Based on this single fit pa- 
rameter P\Qss{Vb), the QT model is therefore a reasonable 
approximation to collisions with unit absorption (or finite 
absorption, by multiplying by an additional absorption 
probability). Moreover, it is easy to evaluate at nonzero 
fields, by simply evaluating Vh for the correct adiabatic 
curve. To a good approximation, the factor Pioss(Vf,) is 
only weakly-dependent on the electric field strength, as 
shown in the inset to Fig. 1. 

This weak dependence can be seen, at least qualita- 
tively, by replacing the actual barrier by an artificial in- 
verse Morse potential, constructed so as to have the same 
curvature at its barrier maximum as the actual potential 
h^L{L+l)/{2^iR^)~Cn/BJ^- The inverse Morse potential 
model can then be solved analytically as a transmission 
proble m [ig, to yield PiQ ss(Vb) = (1 - e-'^''f)/2, where 
/ = \/ L{L + 1) 2(n — 2)/?i . This approximation cor- 
rectly shows that the result is independent of the long- 
range coefficient C„, as well as the reduced mass, but 
that it docs depend on the partial wave / as well as the 
character n of the long-range potential. It also shows 
that, coincidcntally, the transfer function at the barrier 
height is the same for a van der Waals potential n = 6 
and for a dipole potential n = 3, for p-wave collisions. 
Thus a weak dependence of Pioss(Vf,) on electric field is 
perhaps not unexpected. For s-wave collisions, in which 
the long-range dipole potential scales as Ci/R^ [1^, we 
would expect a stronger variation of PiossiVb) with elec- 
tric field. 

Based on these remarks, we turn now to the electric 
field dependence of reactive collisions, making the as- 
sumption that s and y are independent of field. Doing 
so, a numerical calculation readily produces the reaction 
rate constant versus the dipole moment of the collid- 
ing species. Examples are shown in Fig. 2 for identi- 
cal fermions (odd partial waves. Fig. 2a) and identical 
bosons (even partial waves. Fig. 2b). In both cases, the 
overall tendency is for the rates to rise as the field is 
turned on. This rise is, however, more dramatic for iden- 
tical fermions, which arc suppressed in zero-field by the 
van der Waals centrifugal barrier. 

The most striking feature in these figures is the pres- 
ence or absence of resonance-like features as the dipole 
is increased. For weak short-range absorption (e.g., 
y ~ 0.1), these features are pronounced, and fall into 
regular patterns according to the angular momentum L, 
M of the dominant partial wave. For L > these are 
shape resonances behind field-dependent centrifugal bar- 
riers, while for L = these resonances appear as the 
effective long-range potential is systematically deepened 
to include additional bound states [ll|, [l^ . Vice versa, in 
the limit of strong short-range absorption (e.g., y = 1), 
these features are completely washed out. This occurs 
because, at complete unit absorption, the resonant state 
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FIG. 2: (Color online) Dependence of chemical reaction rates 
KiasB on dipole moment d for identical fermionic (a) or bosonic 
(b) KRb molecules. In the case of unit reaction probability, 
2/ = 1, this variation takes a universal form independent of 
details of the short-range physics. For y < 1, non-universal 
resonances appear that reveal more details of the short-range 
interaction. The data from [^ (points) are well-fit by near- 
universal scattering, y = 0.83 



field spectroscopy," [l2|) is a valuable component in accu- 
rately determining the parameters that govern scattering. 
On the other hand, another class of molecules, typified 
by RbCs, are not chemically reactive at low tempera- 
ture [20[, and will yield y = 0, or very nearly so. Here 
the resonances should show up clearly, in cither clastic 
scattering, hypcrfinc-changing collisions, or perhaps in 
three-body losses. Note that the resonances are likely to 
be quite narrow in this case |14{ . 

When the resonances do appear, their positions result 
from phase shifts at short range. We illustrate this in 
Fig. 2b, where resonances are shown for two alternative 
values of s. In this figure the same value of s is used for all 
angular momentum channels {L,M)^ but in fitting real- 
life data this is probably not the case, and each family of 
resonances specified by {L, M) will likely contribute its 
own complex phase shift. Although the overall positions 
of the resonances are not specified, nevertheless their rel- 
ative spacings follow a specific pattern. For example, for 
s-wave resonances the location in field £{n) of the n-th 
resonance takes the form \\% 



decays immediately; it does not survive for even a single 
period of the resonance. 

The presence or absence of these resonances thus 
contains information on the scattering, and in partic- 
ular on the value of y. Consider, for example, KRb 
molecules, which arc chemically reactive at zero temper- 
ature [20, [21| , and should thus possess large absorption 
probabilities y. The points in Fig. 2a show the data 
from the KRb experiment [3 • The best fit value to these 
data yields y ~ 0.83, consistent with large, if not perfect, 
absorption (red line). Interestingly, a fit to the zero- 
field data alone instead yields a value y = 0.4 @. Thus 
the electric- field dependence of reaction rates ( "electric- 



E{n) = E' 



no 



(5) 



where £' is an overall scaling, Urx, is related to the scat- 
tering phase shift at threshold, and no is akin to a short- 
range quantum defect parameter. This approximation 
was derived using semicalssical arguments, which are also 
applicable to the higher-partial wave case [Ij] . The fully 
quantum treatment described here should afford an ac- 
curate, complete parametrization of these features. 
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